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Both non-Abclian gauge fields and minimally interacting masslcss matter fields arc local- 
ized on a domain wall in the five-dimensional spacetime. Field-dependent gauge coupling 
naturally gives a position-dependent coupling to localize non-Abelian gauge fields on the 
domain wall. An economical field content allows us to eliminate a moduli for a instabil- 
ity, and to demonstrate the positivity of the position-dependent coupling in the entire 
', moduli space. Effective Lagrangian similar to the chiral Lagrangian is found with a new 

' feature of different coupling strengths for adjoint and singlet matter that depend on the 

width of the domain wall. 
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Introduction 



Localization of massless gauge fields on a domain wall with 3 + 1 dimensional world volume 

in ■ 

has been a long-standing problem to achieve the dynamical compactification in the brane- 
world scenario [1[ . If the gauge symmetry is unbroken only inside the domain wall and broken 
outside, the gauge field inevitably acquires a mass proportional to the inverse of the width 
of the wall [244] • The localization of gauge fields can be achieved if the gauge theory is in 
the confining phase outside of the domain wall 0, 0] . This requirement can be translated 
into a position-dependent gauge coupling 



The supersymmetric gauge theories in 4 + 1 spacetime dimensions allows a cubic coupling 
0] between adjoint scalar fields S Q and gauge field strengths F^ IN 

C^C a ^ a F^ IN F^ MN , (1.1) 

with coupling constants C a p^. In Ref.Q], a domain wall solution is chosen so that the scalar 
field £ is positive inside and vanishes asymptotically outside of the domain wall. In this 
way, non-Abelian gauge fields have been localized on the domain wall. More recently a 
method to localize non-Abelian gauge fields together with minimally interacting matter 
fields on domain walls has been introduced and a particular model has been presented in 



five-dimensional spacetime [10j] , by gauging the unbroken global symmetry associated to the 



degenerate scalar fields to form the domain wall [ll|]. This mechanism can be viewed as a 
field theoretical realization of D-bran es, and provides a step towards a realistic model of 
branes as soliton solutions of higher dimensional field theories. 




Fig. 1 Position-dependent gauge coupling (right panel) is given as a difference u\ — 02 of 
two separated kinks (left panel). 



In Ref. [Sj, two real scalar fields o\ and 02 were introduced to form the usual kink profile 
in the extra-dimensional coordinate y as shown in the left panel of Fig. [H and the cubic 
coupling in Eq. tjl.ip is chosen as 

£cubic ~ -(<ti - <j 2 )Tt[G M nG m % (1.2) 

where Gmn is the non-Abelian gauge field strength to be localized on the domain wall. 
In this way, one achieved the desired profile of the position-dependent coupling which is 
positive and vanishes asymptotically outside of the domain wall as shown in the right panel 
of Fig(TJ More recently, matter fields have also been localized on the domain wall interacting 



minimally with the localized gauge fields [lOj]. Although the model realized the localization 



mechanism in a simple setup, the width of the position-dependent coupling is a modulus 
and can become negative implying an instability of the gauge kinetic term. 

Moduli of domain walls arise as remaining degrees of freedom of scalar fields constrained 
by gauge symmetry to form the domain wall. This observation prompts us to construct a 
model with smaller number of scalar fields with a different charge assignments in order to 
eliminate the unwanted modulus in the position-dependent coupling Q|. 

The purpose of this paper is to present models in 4 + 1 dimensional spacetime with domain 
wall solutions which allow stable localized gauge fields interacting minimally with localized 
matter fields. For any values of the moduli parameters, we show that the position-dependent 
gauge coupling is positive everywhere in the extra-dimensional coordinate y. In particular, 
its width is no longer a modulus and is fixed by the parameters of the theory. We thus find 
that the gauge kinetic term is stable for all values of moduli, namely for all possible domain 
wall configurations. As localized matter fields, we obtain scalar fields in singlet and adjoint 
representations of the gauge group. They are associated with the broken part of the global 
symmetry, and interact minimally with the localized gauge fields. The effective Lagrangian 



on the domain wall is also worked out, using the method in Ref.|12l|. It resembles the chiral 
Lagrangian associated with the chiral symmetry breaking of QCD. A new feature of the 
effective Lagrangian is that the coupling strength of the adjoint matter fields is larger than 
that of the singlet matter fields. We examine generality of models with a stable kinetic term 
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Table 1 Quantum numbers of fields of the model for the domain wall. 



for gauge fields to identify the class of models with desirable properties. We also discuss 
ways to build more realistic models for a phenomenology of the brane-world scenario. 

The organization of the paper is as follows. Section 2 is devoted to construct a model in 
4+1 spacetime dimensions with domain wall solutions which localizes matter and gauge 
fields. The position-dependent gauge coupling is shown to be positive. In Section 3, the low- 
energy effective Lagrangian is obtained. In Section 4, we examine the generality of models 
with the stable gauge kinetic term. Conclusion and discussion are in Section 5. Some details 
to derive the low-energy effective Lagrangian is given in Appendix A. In Appendix B, we 
summarize geometrical features of the three-flavor model. 



2. A model with localized matter and gauge fields 

In this section we first present a model allowing the domain wall solution with unbroken non- 
Abelian global symmetry [111 ]. As the second step, we introduce non-Abeliangauge fields for 
the unbroken symmetry. As the third step, we consider the cubic coupling [9( in Eq. (|l.ip to 
localize non-Abelian gauge fields : we use expectation values of a singlet scalar field to give 
the position-dependent gauge coupling, whose positivity is demonstrated for any values of 
moduli parameters. A broken part of the global symmetry provides matter fields in singlet 
and adjoint representations of the localized gauge fields 1 1 Of ] . 



2.1. Lagrangian with global symmetry and domain wall solutions 
Let us consider a five-dimensional SU(N) C x U(l)i x £7(1)2 gauge theory and N scalar fields 
Hi (H2) in the fundamental representation with the degenerate mass m (— m), together with 
a singlet scalar field ^3, whose charge assignments are summarized in Tab. [U Therefore 
we obtain global symmetry SU(N)l x SU(N)r x U(1)a- In addition, we introduce adjoint 
and singlet scalars £ and a associated with the gauge group SU(N) C x J7(l)i and U(l)2, 
respectively. We assume the following Lagrangian with the signature (+, — , — , — , — ) 

C = -^(G MN G"») - ±F MN F"» + ^Tr(Av /S ) 2 + ±(d M af 

+ Tr \D M H!\ 2 + Tr\(D M - iA M )H 2 \ 2 + \(d M + iA M )H 3 \ 2 - V , (2.1) 
V = Tr |(E - ml N )Hi\ 2 + Tr |(E - al N )H 2 \ 2 + \aH 3 \ 2 

+ i<7 2 Tr( Cl l^ - HiH{ - H 2 HlY + \(?{c 2 + Tr{H 2 H$) - |F 3 | 2 ) 2 • (2-2) 

The U(N) C = SU(N) C x U(l)i gauge coupling and gauge fields are denoted by g and an 
N x N matrix Wm with M = 0,1,2,3,4. The f7 (1)2 gauge coupling and gauge field are 
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denoted by e and Am- Covariant derivatives and field strengths are defined by 



D M H 1>2 = d M H li2 + iW M H lt2 , D M E = 8 M ^ + % [W M , S] , (2.3) 

and G MN = d M W N - d N W M + i [W M , W N ], F MN = d M A N - d N A M . 

The global symmetry U L G SU(N) L ,U R G SU(N) R , e ia G U(l) A , and the local gauge 
symmetry U c G U(N) C , e 1 ^ G C7" (1)2 ac t on the fields as 

H\ — > e ia U c HiUL , H 2 ^e-^ a+ ^U c H 2 U R , H 3 ^e i0 H 3 , (2.4) 

S -)■ E/cEt^ , cr ^ a . (2.5) 



Let us note that the Lagrangian (|2.2|) can be embedded into the five-dimensional J\f = 
1 supersymmetric gauge theory (with 8 supercharges). This fact allows us to obtain the 
so-called Bogomol'nyi-Prasad-Sommerfield (BPS) domain wall solution as we will show in 
the next subsection. We stress, however, that this particular choice is merely to simplify 
calculations and is not required for our results to hold. 

Without loss of generality, we can assume the mass parameter m to be positive. We also 
assume ci > and c 2 > 0. Then there exist N + 1 discrete vacua with r = 0, 1, ... ,N, where 
scalar fields develop vacuum expectation value (VEV) 

' , #3 = Vc 2 + rci , (2.6) 

(2.7) 

Local gauge symmetry is completely broken and only a subgroup of global symmetry remains 
in these vacua. The breaking patterns in the r = and r = N vacua are 

U(N) C x SU(N) L x SU(N) R x U(l) 2 x U(l) A 

► SU(N) L+C x SU(N) R x U(1) A+C , 

0— th vacuum 

> SU(N) R+C x SU(N) L x l7(l) A -c- 




N—th vacuum 



Let us consider domain wall solutions connecting iV-th (0-th) vacuum at left (right) infinity 
y = —00 (y = 00). Then the (coincident) domain wall solutions preserve the diagonal sub- 
group as the largest global symmetry SU(N)l +r+c , providing the Nambu-Goldstone modes 
associated to the breaking of the global symmetry 

SU(N) L x SU(N) R x U(1) A 

SU(N) L+R+C ' 1 • ] 

We assume that fields depend only on extra-dimensional coordinate y and that all gauge 
fields except W y and A y vanish. 
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Thanks to the special choice of the potential, we can rewrite the energy density as 



£ 



1 



rTr 



D y Z - 9 - (dljy - H X H\ - H 2 H\ 



+ Tr|Z?j /J H" 1 + (S-ml JV )^i| 



+ ^2 (v ~ e2 ( c 2 + ^{H 2 H\ - \H 3 \ 2 )^j + TV \D y H 2 + (S - (a + i^,)l^)ff 2 | 
+ \d y H 3 + (a + ^)^3| 2 + C 2 <9< / CT 



+ a y Tr 



ciS - - mljv) - H 2 Hl(Z - al N ) 



(2.9) 



Thus, we obtain the Bogomol'nyi bound for the total energy (per unit world volume) E 



oo oo 

E= j dy£>T = ci j dy d y Tr(T,) = Nma , 



(2.10) 



where T is the tension of the domain wall. This bound is saturated when the following BPS 
equations are satisfied 



d y H x + (S + iW y - ml N )Hx = , 
d y H 2 + (s + iWj, - (a + iA y )l N ^jH 2 = , 
d y H 3 + {a + iA y )H 3 = 0. 



8 y a = e 2 (c 2 + Ti{H 2 H\) - \H 3 



(2.11) 
(2.12) 
(2.13) 

(2.14) 
(2.15) 



To use the moduli-matrix formalism [131-4151] . we introduce S(y) G GL(N,C) and ijj(y) G 



S + iW y = S- x d y S, a + iA y = -dyijj . 



(2.16) 



Then the matter part (|2.1ip - (|2,13p can be solved by 



H 2 = e^S~ l H% , 



H 3 



(2.17) 
(2.18) 

'3 , (2-19) 

with complex constant N x N moduli matrices H® , H® and a moduli constant H® , which 
describe moduli of the solution. The rest of the BPS equations (|2.14[) and (|2.15p turn into 
the master equations for gauge-invariant Hermitian fields O = SS* and rj = ~Re(ijj). 



dyidyQn^ 1 ) = ^ 2 (ci1tv - {e 2my H^H^ + e^H^H^)^ 1 
Itfr, = e 2 (c 2 + ^TriH^n- 1 ) - e~* \H%\ 2 



(2.20) 
(2.21) 



Moduli matrices related by the following ^/-transformations give identical physical fields 

(S, V, Hi Hi Hi) -> (VS, i; + v, VHl VH° 2 e-^,H° 3 e l n , (2.22) 

where V G GL(N, C) and v G C. The equivalence class quotiented by this ^-transformation 
defines the moduli space of domain walls. We can use this freedom to choose the form of the 
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moduli matrices 



B% = yfcl Nl H° 3 = ^. (2.23) 

Let us also decompose H® as 

Hi = yfceW , (2.24) 

where <j> is a Hermitian N x TV matrix and U is a unitary N x N matrix. With this choice, 
the master equations (|2.2U|) and (|2.21|) become 

dyidySin- 1 ) = Un - non- 1 ) , (2.25) 



±d 2 V = e 2 (c 2 + c ie "Tr( e 2 ^- 1 ) - e^pa) , (2.26) 

where ft = e^l^ + e^e". 

No analytic solution of this system of the differential equations is known in general. How- 
ever, one can study essential features of solutions, if one takes the strong gauge coupling 
limit <? 2 ,e 2 — > oo. Eqs. (|2.25p and (|2.26|) reduce to a system of algebraic equations in this 
limit: 

U = e 2my l N + e^e* , (2.27) 
C2 = _ Cie r 'Tr(e 2 ^fi- 1 ) . (2.28) 

It turns out that the effective theory describing massless excitations localized on the back- 
ground solution of the equations of this system precisely coincides (at least in the lowest 
order of approximation) with the one obtained from (|2.25p and (|2.26p (see the detailed 
discussion in Ref. [lOj). It is therefore sufficient just to study solutions of ([2.27P and (|2.28p . 



Without loss of generality we can assume that the moduli <j) is diagonalizable 

= mP- 1 diag(y 1 ,...,y JV )P. (2.29) 
Then (|2.28p reduces to a polynomial equation of order N + 1 for x := e~ v 

x = l + ^Y , e t = e 2m ^\ (2.30) 

C2 f-f 1 + etx 

i=i 

If this equation is solved, one can supply its solution into Eq. fl2.27j> to obtain Q. 

In the simplest case, where all walls are coincident <j) = myolw, we can solve equation 
(f2~30l) explicitly (e := e 2m ( y - y °^) to find 

e~ v = ^ (e - 1 + V(l-eo) 2 + 4(l + iV Cl / C2 ) e o) , (2.31) 
n = (e 2my + e 2myo e r >)l N . (2.32) 
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Fig. 2 Profiles of e v and a in the coincident case. The parameters of the plot are given 
above the picture. Positions of all domain walls are centered at the origin. 

Physical fields can be expressed in terms of f2 and a as 

H, = ycl . ljV , == , (2-33) 

yj\ -\- e -2m(y-y )+r] 

H 2 = , , (2-34) 

#3 = ^e""/ 2 , (2.35) 

£ = -$,lnn, (2.36) 

a = d y rj , (2.37) 
W„ = A y = , (2.38) 

where we fixed the gauge such that 5 = fi 1 / 2 and Im(^) = 0. 

This set of solutions is not invariant under the symmetry transformations (j2.4j) in general. 
For a choice of U = In, however, the solutions (I2.33p - (|2.38p are invariant under the action 
of the diagonal global symmetry SU(N)l+r+ c - We show the y-dependence of e~ v and of 
g = dyTj/2 for the coincident case in FigJ^l 

For more general case such as non-coincident walls, the dependence of e _T? and a on y is 
more complicated. Furthermore, the equation (|2.30p cannot be solved in a closed form in 
general, except for first few values of N. Thus, one has to use numerical techniques. In Fig. 
[31 we present an example of five non-coincident walls. 

2.2. Localization of non-Abelian gauge fields 

In order to obtain massless gauge fields localized on the domain wall we need to introduce 
a new gauge symmetry which is not broken in the bulk. As we have seen in the previous 
subsections, the coincident domain wall solutions fl2.33jl - fl2.38D do not break a large part of 
the global symmetry. Let us then gauge SU(N)l + ji = SU(N)v and denote new gauge fields 
as Vm- Then the fields H\ and H2 are in the bi-fundamental representation of SU(N) C x 
SU(N)v and the covariant derivatives fl2.3h are modified to 

£>m#i,2 = d M H 1>2 + iW M H lj2 - iH lt2 V M ■ (2.39) 
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Fig. 3 Profiles of e ^ and a in the non-coincident case. The parameters of the plot are 
given above the picture. Positions of domain walls are y\ = —10, yi = 4, y% = 8, 2/4 = 12 and 
2/5 = 16. 
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Table 2 Quantum numbers of the gauged model. 



Quantum numbers of the gauged model are summarized in Tab. [2] 
We introduce field-dependent gauge coupling for Vm as 

1 



Act , (2.40) 



where we assume that A is a positive constant 

A > 0. (2.41) 

If we denote the field strength for the new gauge fields as Gmn the Lagrangian for the 
gauged model is given by 



C = I — Tr 

23 V) 



GmnG 



AIN 



(2.42) 



where C is the same as in (|2.ip except that the covariant derivatives (|2.3f) are replaced 
by (|2.39p . The choice of the field-dependent coupling (|2,40p is inspired by supersymmetry. 
As discussed in Ref. 0, a term that is linear in adjoint scalars, appearing in front of the 
kinetic term for gauge field, naturally arises in five-dimensional J\f = 1 supersymmetric gauge 
theories & In tat, our model can be embedded into M = 1 supersymmetry. 

It is not hard to see that the solution f|2.33j) - (|2.38[) in the ungauged model is equally valid 
in the gauged model. If we write down the equation of motion for the new gauge field Vm 
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we have 

8 M G MN = J N , (2.43) 

where Jm stands for the current of Vm- Since the solution preserves SU(N)v, the current 
vanishes for the domain wall solution (|2.33|) - (j2.38p . and Vm = is a valid solution to the 
equation of motion (|2.43p . Then the other equations of motion of the gauged model reduce 
to those of the ungauged model because of Vm = 0. Therefore, we see that f)2.33j) - f|2.38j) in 
addition to the condition Vm = solves the whole set of equations of motion of the gauged 
model. 

2.3. Positivity of position- dependent gauge coupling 

We wish to show that the field-dependent gauge coupling (|2.40p assures the positive defi- 
niteness of the position-dependent gauge coupling for any configurations of the domain- wall. 
Since we do not need the effective theory in full, we reserve to derive the rest of the effective 
Lagrangian for the next section. 

For the moment it is sufficient to know, that the field-dependent gauge coupling l/g 2 (a) 
is given by its value in the background solution 

\d y n = -Xdylnx, (2.44) 



background 



~g 2 (y) 



where x = e ^ > is a solution to (12.301) . Differentiating (I2.30P we find 



1 d„x = - — - — e% , 9 (2m + -d y x ) . (2.45) 
c 2 (1 + eix) 2 \ x y J 



v 

x y 



This leads to the formula 

WU) = ^ § (1 + ^) 2 / i 1 + ~2 § (l + e^) ' (2 ' 46) 

which is indeed positive in the whole range of y-coordinate. 

Integrating l/g 2 (y) over the extra-dimensional coordinate y we obtain the effective gauge 
coupling in 3 + 1-dimensional world volume 

oo 

^2=^j dyd yV = X[ V (oo)-7 ] (-oo)]. (2.47) 

— oo 

The asymptotic values of r\ are found from (|2.3ip as 

77(00) = 0, r?(-oo) = - In M + V (2.48) 

The easiest way to see these asymptotic values is to note that in (|2.32p we can take the 
limits of (|2.3Up to obtain 

N 



since g is finite at both infinities 0. 



We can just look at (|2 . 19[) with = ^foi and recall (|2.6p to obtain the same result. 
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Fig. 4 Profile of ?7-kink is shown in the left panel for the coincident case. In the right 
panel, plots of Tr[£] (green dashed curve), Tr[E] — a (red solid curve), and a (blue dotted 
curve) are shown. 



Thus, the effective gauge coupling is given as 

-4 = Alnf 1 + iV— ] . (2.50) 
9 \ c 2 J 

It is interesting to observe that the effective gauge coupling is proportional to the width of 
the domain wall ln(l + NC1/C2), which is not a modulus, but is fixed by parameters of the 



theory. This feature is in sharp contrast to that in Ref. [10(j where the effective gauge coupling 
constant is proportional to the domain wall width, which is a modulus undetermined by the 
theory. We have now confirmed the stability of the gauge kinetic term (by choosing the 
parameters of the theory as A,m, c±, oi > 0). 

Eq. (|2.47p shows that the effective four-dimensional gauge coupling constant g is deter- 
mined only by the boundary conditions at infinity. It can be interpreted as a kind of 
topological charge of the 77-kink, whose profile is shown in Fig. HI 

Note that we have considered only the BPS solutions so far, and confirmed the stability of 
the model. The anti-BPS solutions are also stable, since they can are obtained by the parity 
transformation x — > — x. 



3. Effective Lagrangian 

In this section we calculate the low-energy effective Lagrangian on the background domain 
wall solution in the moduli approximation 16], where the moduli are promoted to fields on 
the world- volume with coordinates and are assumed to depend only weakly on x^. As the 
background we consider the coincident walls 4>b g = myolv- Our model have two moduli: a 
Hermitian N x N matrix <j) and a unitary N x N matrix U. We study in two steps: taking 
only U as moduli fields, and then both <f> as well as U. 



3.1. Effective Lagrangian for coincident walls 

Here we ignore the moduli fields <p, and take only U moduli fields, that describe the Nambu- 
Goldstone modes associated with the symmetry breaking in Eq. (|2.8p . When <j) = myoljv, 
the solution (|2.17p - (|2.19|) with (|2.23|) and (|2.24|) become in the strong gauge coupling limit 
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g — > oo and e — > oo as 



H 2 = ^e^e^^n-^U^x) , 



c 2 e 



-rj/2 



(3.1) 

(3.2) 
(3.3) 



where O and are denned in Eqs. (|2.31|) and (|2.32|) with the replacement yo ~~ * 2/o( a;/i ) an d 
£/ — >■ Ui^x 11 ). We plug these into the Lagrangian (I2.42H . where the gauge fields and A M 
are no longer dynamical and should be eliminated as auxiliary fields. After integrating over 
the extra-dimensional coordinate and taking up to quadratic terms in the derivatives, we 
obtain the low energy effective Lagrangian (the detailed calculation is given in Appendix lAl) 



Left — ~ — 



+ 



2m 
Nmc\ 



(q + lyUrlpuZflTPU] + — Tr[UV^]Tr[UV^) 



d^yod^ya - ln(l + — ) 



C2 



where 



and 



v^ = d^u + i[v^,u], 



1 , c 2 c 2 
a = — 

2 N Cl Na 



1 + 



In 1 + 



Nci J \ c 2 

Let us define the decay constants f n for the adjoint field and f v for the singlet field 



(3.4) 

(3.5) 
(3.6) 



ci(a + 1) 
2m 



fv 



Nm 



(3.7) 



Then the canonically normalized adjoint and singlet fields can be defined as tt and rj 
respectively 



— £>»7T 



Jr] ^ ' 

The effective Lagrangian Eq. (|3.4p can be rewritten as 



(3.8) 
(3.9) 

(3.10) 



A new feature of (I3.10|) compared with our previous work [10| is that the coupling strength 
f n of the adjoint field is larger Q than /„ of the singlet field by a factor \/a + 1. 



t Difference of factors of 2 and N is due to a convention of SU(N) matrix normalization 
Tr(T I T J )=S IJ /2. 
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3.2. Full effective Lagrangian 

As the final step, we consider the general case of both U and d> as moduli fields. Then scalar 
fields are given in terms of moduli fields U(x fl ) and cj){x^) as^j 

Hi = ^e my n~ 1/2 , (3.11) 

H 2 = ^/cle^ir^Vt/t , (3.12) 



c 2 e 



-rill 



n- 1 



e -v e -W 



1 + ^Tr 



(3.13) 
(3.14) 

\1 N + e 2ray e -2<t> e -ri I 

To obtain the effective Lagrangian, we need to repeat the same procedure as in the previous 
subsection, where the covariant derivatives acting on functions of matrices require more care 
(see e.g. Appendix B in Ref. lOJ]) and causes difficulty when deriving the closed form of 
the effective Lagrangian. However, we have a convenient parameter to expand the effective 
Lagrangian, the ratio c\/c 2 whose logarithm has a physical meaning as the width of the 
domain wall (see Eq. (|2,5U|) ). As given in Appendix the effective Lagrangian up to the 
order of c\ / c 2 is given as 

A* = 4? + T} 1] + T$> + T^i + % + % + Cl 0({ Cl /c 2 f 



(3.16) 



where 



cff 2m 



<b C0 . Sh{C f ) 7 1 In ( 1 + tanh(£ ^' ] (DM, 



+ 



sinh(£^) \ 1 - tanh^) 

c°Bh(£,) - l ^ / l + U^j 
Ccj, smh(£^) Vl-tanh(£^)/ 

2 ^ tanh^) 



,l-tanh(£ )7 v 'J' (3 ' 17) 

with a Lie derivative with respect to A, Ca(B) = [A, B\. Interestingly, this is the same 
effective Lagrangian we obtained previously [10|. The rest of terms are of order of c\/c 2 , 
given by 



T, 



(i) 



u 



I6c 2 m 



r 



(i) 



,2 
c l 



8c 2 m 



(i) 



8c 2 m 



: Ti[v^UFu(d x ,C^(WV^U)e 2 
U^UFrnUd^jC^CD^e**} Tr 



Tr 



-X(f> 



x=0 



■Tr 



x=0 



-Tr 



% 



lGc 2 m 



F(d x )Tr 



Tr 



I6c 2 m 



F(O x )Tr 



Tr 



Tr 



x=0 



(3.18) 
(3.19) 
(3.20) 
(3.21) 
(3.22) 



x=0 



* The center of mass position of walls yo(x fl ) is a free field as in Eq. (|3.10p and is suppressed here. 
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and 



DO 



Fu(x C)= I dy cosh(£^) 

cosh 2 (y — x) cosh(y) cosh(y — C^) ' 



-oo 
oo 



cosh(y) coshfu — £&) — sinh(y) sinh(y — C&) — 1 , 

F mix (x,C^= dy — ^- — ^ ^ , (3.24) 

£^ cosh(y) cosh(y — C^) cosh (y — x) 



-oo 
oo 



-oo 
oo 



, cosh(y) cosh(y + - sinh(y) sinh(y + - 1 

F (f> {x,C 4> ) = / — -j- f , (3.25) 

LA cosh(y) cosh(y + L^) cosh [y — x) 



-oo 



1 

cosh 2 (y) cosh^(y — x) 



F(x)= I dy—^- -2—— -. (3.26) 



All the above integrals can be obtained in closed forms. 

4. More general models of stable position-dependent coupling 

In this section, we wish to show that there are more models with the stable position- 
dependent coupling. We will illustrate the point by extending the model to include more 
fields and more gauge symmetry. 

The position-dependent gauge coupling comes from the cubic coupling between a singlet 
scalar field and field strengths of non-Abelian gauge fields in 4 + 1 dimensions 

Aubic = CV 4 )Tr [GmnG mn ] , (4.1) 

where the function C(<7j) of singlet scalar fields Oi should be linear, if it is to be embeddable 
into a super symmetric gauge theory in 4 + 1 dimensions 

C(ai) = J2^h 7«£K, (4.2) 

i 

where 7, are constant coefficients. 

Usually each domain wall has one complex moduli: a position and a phase. For exam- 
ple, both the massive CP 2 model and the massive CP 1 x CP 1 model have two free domain 
walls, corresponding to the two complex moduli. Although two free domain walls can pro- 
duce the desired profile of position-dependent gauge coupling by an appropriate choice of 
parameters in Eq. ()4.2p . they provide the undesired modulus for the width of the profile. 
To avoid this problem, we are led to consider models with a single complex moduli. The 
simplest one of such models is the three flavor model in Ref.Q], where three scalar fields are 
constrained by the two Abelian gauge symmetry U(l) x U(l). Geometry of this three- flavor 
model is examined in Appendix [Bj Our model in this paper is an extension of this model 
to non-Abelian gauge group: U(l) X U(l) — > U(N) x U(l). The next-simplest possibility is 
to consider four scalars constrained by three Abelian gauge symmetry U(l) x U{1) x C7(l) , 
which we call four-flavor models. We give quantum numbers of fields of a typical four-flavor 
model in TabJH In the limit of strong gauge couplings, the gauge theory becomes a nonlinear 
sigma model whose target space is given by an intersection of three conditions as 

(c 2 x cp 1 ) n (m 2 x cp 1 ) n (c 2 x cp 1 ) ~ CP 1 . (4.3) 
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CCD] 


^ \ J* 


U(l) 3 


mass 




1 








TYl\ 


Ho 


1 


-i 





ui2 


H 3 





1 


1 


m 3 


H 4 








1 


1714 


<y\ 














02 














0-3 















Table 3 Quantum numbers of the U(l)i x U(l)2 x £^(1)3 four-flavor model. 





W 


\H 2 \ 


\H 3 \ 


W 


0"1 


&2 


0-3 


(1) 





Vci 


V /C l+2 




m2 + m 3 — 1714 


m 3 — m A 


7714 


(2) 









V C 3-2 


mi 


m.3 — 7714 


7714 


(3) 









\/C3 


mi 


mi — m 2 


7TI4 


(4) 


V C l+2-3 


\/ C 3-2 


\/C3 





mi 


mi — m 2 


— mi + 7712 + ^3 



Table 4 VEVs of candidate vacua: We use abbreviations like C3_i_2 = c 3 — ci — C2. 



The vacuum condition is given by 

\Hi\ 2 + \H 2 \ 2 = ci, -\H 2 \ 2 + \H 3 \ 2 = c 2 , \H 3 \ 2 + \H 4 \ 2 = c 3 , (4.4) 

H x {gi - mi) = 0, ff 2 (o"i - cr 2 - m 2 ) = 0, (4.5) 

#3(0-2 + 03 - m 3 ) = 0, #4(03 - m 4 ) = 0, (4.6) 

where c% is the Fayet-Iliopoulos parameter of U(l)i and m a is the mass for Hi. All possible 
solutions to these equations are shown in Tab. [H There are four solutions but only two of 
them are valid solutions for any choice of real parameters of Cj. When we choose ci > 0, 
C2 > and c 3 > ci + c 2 , we are left with the vacua (1) and (2) in Tab. [U 

The moduli matrix formalism 0, [TtJ is powerful enough to give generic solutions of the 
BPS equations of this nonlinear sigma model. Especially, we are interested in the kink profiles 
of <7i , o"2 and 03 . They can be expressed as derivatives of real functions r\i 

o-i = -d y rn, (i = 1,2,3), (4.7) 
where the real functions r\i are determined by the following algebraic conditions 

e -r?i+2mia _|_ e -r) 1 +r}2+2m 2 y-2a _ £ ^ /^_g\ 
_ e -r]i+V2+2m 2 y-2a _j_ e -77 2 -r?3+2m 3 y _ ^ /^g\ 
e -r? 2 -773+2,^33/ + e -%+2m 4 y = ^ ( 410 ) 

with a real constant a. The parameter a is (the real part of) the unique modulus of the 
solution, corresponding to the position of the domain wall. As we expected, we find only a 
single modulus. The width of the domain wall is not a modulus, but fixed by the theory. 

Since we want a configuration that a% — > at both spacial infinities, we can choose m 3 = 
777,4 = 0. Then 02 = 03 — > at y = ±00. Several numerical solutions are displayed in Fig. 
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Fig. 5 The kink profiles of o\ (red solid line), 02 (green dashed line) and 03 (blue dotted 
line) for two different sets of model parameters. 



From Fig. [5l we clearly see that 02 > and 03 < for all the values of y. As a position- 
dependent gauge coupling, we can choose a two-parameter family of desirable models 

Ccs = ~ (7202 - 73^3) Tr [g mn G mn ] , (4.11) 

where 72,3 can be any non-negative real numbers. This class of models can be easily made to 
localize non-Abelian gauge fields and minimally interacting matter fields by extending two of 
the U(l) factor groups to (possibly different) U(N) gauge groups with N scalar fields in the 
fundamental representations, similarly to our model in previous sections. A new interesting 
feature of the four-flavor model (and its non-Abelian extensions) is that two possible profile 
of singlet fields u\ , 02 can provide a different profile for different non-Abelian gauge groups 
such as SU(3), SU(2) and U(l) with their associated matter fields. 

5. Conclusion and discussion 

We have successfully stabilized gauge fields and matter fields that are localized on domain 
walls. The low-energy effective Lagrangian on the domain wall has been worked out, where 
the adjoint matter fields is found to couple more strongly than singlet matter fields. We also 
explored possible generalization of our stabilization mechanism by including more fields and 
more gauge symmetries and found a class of more generic models with an added flexibility 
for model building with different localization profile for different gauge groups. 

To build realistic models of brane-world with our scenario of localized gauge fields and 
matter fields, we should address several questions. Perhaps the most important question 
is to obtain (massless) matter fields in representations like the fundamental rather than 
the adjoint of localized gauge fields. One immediate possibility is to use the localization 
mechanism of fermions in a kink background. It has been found that zero modes of such 



fermions are localized in such a way to give automatically chiral fermions 18(] . We will pursue 
this direction and associated anomaly questions further. 

Secondly, we should devise a way to give small masses to our matter fields in order to 
do phenomenology. Since some of our matter fields are the Nambu-Goldstone modes of a 
broken global symmetry, we need to consider an explicit breaking of such global symmetry. 
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Thirdly, another question is to study possibility of supersymmetric model of gauge field 
localization. We need to settle the issues of possible new moduli in that case. Moreover, we 
should examine the mechanism of supersymmetry breaking. 

Finally it is an interesting possibility for model building to localize gauge fields of different 
gauge groups with different profiles. This situation is often proposed in recent brane-world 
phenomenology, for instance in Ref. 19 ] . 
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A. Derivation of the effective Lagrangian 

In this appendix, we derive the effective Lagrangians (|3.4f) and A3. 161) in the strong coupling 
limit. Since we are interested in the low energy effective Lagrangian, we focus on up the 
quadratic terms in derivatives, which we call £( 2 \ 



A.l. Eliminating gauge fields 

Here we eliminate gauge fields W^,A^ to obtain a nonlinear sigma model. Starting with £ 
in (j2.42p in the strong gauge coupling limit g — > oo and e — > oo, is given as 

=Tr| J D^ 1 | 2 + Tr|( J D At -a M )F 2 | 2 + |(^ + ^ At )F 3 | 2 , (Al) 

where the covariant derivatives are given as: 

D fl Hi,2 = Au#i,2 + = d„H 1)2 - iH^ + iW„H lj2 . (A2) 

Here we singled-out covariant derivatives containing fields associated with the gauged 
part of the flavor symmetry SU(N)l+r- Notice that Wn and are no longer dynamical, 
but they are merely auxiliary fields to be eliminated. Moreover, the constraints are satisfied 
by the scalar fields (|2.33|) -( i2.35p which depend on the moduli fields 4>{x^) and U(x fl ) 

H 1 hI + H 2 hI = c 1 1 n , (A3) 
\H 3 \ 2 -Tr(H 2 Hl) = c 2 . (A4) 



The equation of motion for Wn gives 



I 

2c7 



HaD^Hl - t)^H a Hl] + -A^H 2 H j 2 = W^ + -A^H 2 Hl , (A5) 

Cl Cl 



1G 



where the sum over the index a = 1, 2 is implied (Einstein summation convention). Plugging 
this back into (|Alj) we obtain 

+ A^(\H 3 \ 2 + Ti(H 2 Hl) - ±-Tr[(H 2 Hl) 2 ]) + d^H^Hs . (A6) 

Next step is to eliminate the auxiliary fields A^. By using the following identities derived 
from the constraint ( 1A3D 

(A7) 
(A8) 



Tv[(H 2 Hl) 2 ] = Cl Tt(H 2 Hl) - TviH^l^Hl) , 
-2Tt{W»H 2 hI) = iTr(H 2 D^Hl - Z^feflJ) 

+ -eabTrliHaD^Hl - D„H a Ht)H b Hl] , 

and by solving equations of motion for A^, we obtain 

i H^H\ - d^H 3 Hl + j-e ab Tv[(H a D^Hl - b^H a H\)H b Hl] 



A,. 



\H 3 \ 2 + ^(HvHlH^ 



(A9) 



By using the following identities 



Tr 



2ci 



where U a b ■= iAn h and V^H ab = d^H ab + i[V^, H ab ], and 

eabTrkHaDpHt - D^H a Hl)H b Hl] = Tx[u\ 2 V ^ - V^U^U 
we obtain a simpler expression for A^ 

H z d„Hl - d^H 3 Hl + j-Tr[u\ 2 V^U 12 - V^u\ 2 )U l2 



12 



-4, 



\H 3 \ 2 + j-Tr[H\ 2 H 



(A10) 



(All) 



(A12) 



12 



Using (IA12j) . we can rewrite the effective Lagrangian as an integral of a nonlinear sigma 
model over y 



£ c ff 



/ 



dy 



1 

2c~i 



Tr 



V,U ab V»-H\ b 



+ d„H 3 d»H. 



t 



^(|F 3 | 2 + iTr(^ 12 )) 



(A13) 



t4.£. Effective Lagrangian for U and yo 

Let us now calculate the effective Lagrangian including only fluctuations U and yo around 
the coincident domain wall solutions (|3.1|) - (|3.3p with (|2.3ip and (|2.32p . The composite fields 
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T~L a b are given as (eo = e 

^11 = 01 ° 1 N , (A14) 
1 + e e ^ 

e V2 -„/2 

«i2 = ci i e ° (A15) 
1 + e e 

^22 = ci 1 - Ijv- (A16) 
1 + e e 

After some algebra, covariant derivatives of these can be rewritten as 

VJlii = -^e-V^yoliv , (A17) 
T>^Hi2 = ^e~Vsinh(m(y - y Q ) - i]/2)d^yQ U ] 



ci 



+ y w ; * (A18) 

2 cosh(m(y -y ) -r]/2) 

V^n 2 2 = ^e-V^yoljv , (A19) 



where 



.l + -^(l + e oe ^) 2 /eo ■ (A20) 
m V Wei / 



Putting this into Eq. (|A"T2|) we obtain: 



^ = "2^^^^ " • (A21) 



Substituting this result into (IA13p . we obtain 

oo oo 

I dy-^-j-e- T iTT[V^V^U]+mc 2 J dye~ r 'ad^y Q 



Nm 



-oo 

oo 



+ / dV ^j^e-^iUV^MUV^] . (A22) 



-oc 



Finally integration over y is most easily done by substituting x = e ri and using the identity: 



x + (x-l)(l-^-(x-l) 



(A23) 



leading to Eq. (|3.4p (including the kinetic term for the gauge fields). 
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A. 3. Effective Lagrangian for all terms: <f>, U and yo 

Next we derive the effective Lagrangian including all fluctuations. Plugging the solution 
(I3.11|] - (I3.13|) into the Lagrangian (]A13p . we obtain 



-Ccff 



dy 



Tu + Tmix + %> 



+ 



16c 2 



dy e*(l-—)Tr 



-oo 
oo 



16c 2 



mJ 



a 



1 



dye" 1 Tr 

m 



cosh 2 (y) 
1 

cosh 2 (y) 



V„U ] U 



Tr 



1 



Tr 



cosh 2 (y) 
1 



cosh^(y) 



(A24) 



where 



4 [ cosh 



«) 



\ (") p-y \ 



(A25) 



77, 



ci ^ ("l) r 



n=2 



(«) 



1 \ (n) 1 



cosh(y)7 cosh(y) 



+ ( tanh(y) I tanh(y) 



(A26) 



and 



4 1 — a/m 4 



1 V^V^ 
cosh (y) 



n=3 

4- I tanh(y) I tanh(y) 



vmv»<i> 



l 



CO 



1 



Acosh(y)/ cosh(y) 



(A27) 



y = (my- V /2)l N -(f>. (A28) 

The only work which remains to be done is to perform the integration. This, however, turns 
out to be difficult. The main source of difficulty lies in the fact the we cannot solve ()3. 15[) 
explicitly. We therefore use some approximation technique, such as the Taylor expansion. 
We hnd it convenient to use as an expansion parameter c\/c2 which determines the width 
of the coincident wall. In the lowest order of approximation we see that Eq. ()3.15p reduces 
to e _?? = 1 and we obtain the effective Lagrangian in Eq. (j3.16[) . 



A. 4- Complexity of the Effective Lagrangian 

The formula (|3.16p illustrates the complexity of the interactions between moduli fields <j) 
and U in the general case. Let us offer some explanation for this complexity. The structure, 
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Hi 





-1 





cr\ 











02 












Table Bl Quantum numbers of the U{\)i x U{1) 2 three-flavor model. 



which is new in (I3.16P compared to (I3.17P is of the general form 

F(^)Tr[e^M(y)]Tr[e-^iV(y)] , (A29) 

x=0 

where M(y) and N(y) are some matrix-valued functions, containing either derivatives of cf> 
or derivatives of U . If we assume that <fi = mP~ 1 diag(yi, . . . , dn)P is diagonalizable we can 
rewrite the above as 

N 

£ F(m( yi - y J ))(PM(y)p- 1 ) ll (PiV(y)p- 1 ) ij . (A30) 

This form suggests that interaction (here represented by function F) depends on the relative 
size of fluctuation of each pairs of walls. Indeed, notice that if two walls have the same 
position yi = yj i ^ j (meaning that the expectation values of the fluctuations is the same), 
the above form reduces to 

F(0)Tr(M)Tr(iV) , (A31) 

which is in a sense trivial, since we already encountered this kind of terms in (|3.4p . Thus, 
the new kind of complexity in our result (|3.16|) can be understood as a manifestation of the 
fact, that the interaction does not depend only on various moments of the fluctuation as in 
(I3.17D but also on their relative size. 

B. Geometry of the three-flavor model 

Quantum numbers of three-flavor model are shown in Tab. IBll Since we are interested in 
the domain wall solutions, it is enough to consider the strong gauge coupling limit where 
gauge theories become non-linear sigma models (NLSM), whose target space is defined as 
an intersection of two spaces as 

(C x CP 1 ) n (C x M 2 ) ~ CP 1 . (Bl) 

Here H 2 stands for the two dimensional hyperbolic plane. In the above expression, The CP 1 
and M 2 are defined by 

CP 1 = {(Hi,H 2 ) | |iTi| 2 + |# 2 | 2 = ci, Hi t2 eC}/U(l)i, (B2) 

H 2 = {(PT 2 ,Pf 3 ) | \H 2 \ 2 - \H 3 \ 2 = -c 2 , %€C}/f/(l) 2 . (B3) 

This space is isomorphic to CP 1 but it is a squashed sphere. The metric can be read from 
the Kahler potential 

K = e~ v '\Hi\ 2 + e- v ^ v *\H 2 \ 2 + e v *\H 3 \ 2 + c x V x - c 2 V 2 . (B4) 
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One can eliminate the real superfields V\ and V2, and find the following expression of the 
Kahler potential with respect to the gauge invariant inhomogeneous coordinate cp as 

K = dfj + logtM^ + r^-Alog/], (B5) 
1 



/ = ^A-M 2 + y4M2(l + A) + (A-M2) 2 J, (B6) 

H 2 H 3 c 2 . , 

99 = — — , A=— . B7 
C1-H1 ci 

Note that this manifold is singular at 99 = when C2 = 0. This can be understood in two 
different ways. The first one is to realize the fact that ^(1)2 gauge symmetry is restored 
at that point. This is because = \H%\ = holds there. Namely, the additional massless 
degrees of freedom should be taken into account. The second way is more straightforward. 
Let us calculate the scalar curvature in the vicinity of (p = 0. To this end, first we change 
the coordinate by <p = e** tan ^| (0 < 6 < 7r, < < 2tt), then we obtain 

fl= 8(8 + 9A) 1 + 
9 c 2 

From this it is clear that the scalar curvature at (p = becomes infinity when c 2 = 0. 
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